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Abstract 

The existence of radial solutions of a nonlinear Dirichlet problem in 
a ball is translated to the language of Mechanics, i.e. to requirements 
on the time of motion of a particle in an external potential and under 
the action of a viscosity force. This approach reproduces existing 
theorems and, in principle, provides a method for the analysis of the 
general case. Examples of new theorems are given, which prove the 
usefulness of this qualitative method. 
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1 Introduction 



In the present paper, we consider the following nonlinear Dirichlet problem 



where / is a differentiable function and Q is the ball of radius R in $l D . 
We look for conditions guaranteing the existence of spherically symmetric 
solutions to (1-2). 

The above mentioned problem has been extensively studied in the last 
years (see, for example, [1-5] and references therein). In this paper, our pur- 
pose is to develop a very simple picture, based on Mechanics, for the analysis 
of the existence of solutions to (1-2). This qualitative picture reproduces 
the existing results and, in principle, provides a frame for the analysis of 
the radial solutions to (1 - 2) in presence of an arbitrary nonlinear function 
/. Examples of new theorems are given, which show the usefulness of the 
method. 

To our knowledge, the analogy of the radial equation (1) with the Newto- 
nian law of motion of a particle was first used by Coleman [6] to obtain the 
approximate form of the solution connecting false and true vacua in scalar 
field theories. This solution enters the semiclassical expresion for the decay 
probability of the false vacuum state. Application of this analogy to the anal- 
ysis of the existence of solitary waves in nonlinear one-dimensional media has 
proven to be very useful too [7]. 

The plan of the paper is as follows. In the next Section, the problem about 
the existence of solutions to (1-2) is translated to the language of Mechan- 
ics. Two limiting solvable cases, the one-dimensional problem and the linear 
equation, are considered and a few general results are given. Let us stress 
that the function f(u) is interpreted as the derivative of a potential, thus the 
linear equation describes the motion in a quadratic potential. Section 3 deals 
with potentials having a well around u — 0. The most interesting examples 
studied in this Section are, in our opinion, the potentials with barriers. In 
Section 4, we study the motion in a potential with a hill around u — 0. In 
Section 5, we consider singular (finite and infinite) potentials. Concluding 
remarks are given at the end of the paper. 



Au + f(u) = in Q, 
u = on <9f2, 



(1) 
(2) 
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2 The Analogy with Mechanics 



We start by considering the spherically symmetric version of Problem (1-2) 



d 2 u D — 1 du 



dr2 +— Tr +/(«) = 0. P) 

fill 

^(0) = 0, u(i2)=0. (4) 

Written in this form, the analogy with Mechanics is evident. Equation 
3 is nothing, but the Newton law for a particle of unit mass moving in a 
potential^ which is the antiderivative of /, f(u) = dV/du, and under the 
action of a viscosity force inversely proportional to time. The particle should 
start with zero velocity from a position u(0) and arrive to u = in a time R 
(Fig. la). 

We have drawn in Fig. lb a generic positive solution to (3-4) for a given 
V. In general, the particle will realize damped oscillations around the point 
u — (Fig. 2). Let T n (u(0)) be the time the particle spends to reach the 
point u = n times starting from u(0). Thus, the existence of a solution to 
(3-4) may be formulated in the following terms: 

"In the potential V, there exists an u(0) and a positive integer, n, such 
that T n (u(0)) = R" 

The interesting point is that in many cases we may perform simple es- 
timates, based on physical principles, of the dependence T n vs u(0) and, 
consequently, we may give criteria for the existence of solutions to (3-4). 

Let us first study two limiting cases in which equation (3) may be solved 
exactly. They will be very useful in the analysis below. 



2.1 The one-dimensional (D = 1) problem 

The D = 1 case is characterized by the absence of friction. Thus, the energy 
E = {l/2){du/dr) 2 + V(u) is conserved, dE/dr = 0, and the dependence r{u) 
may be expressed in the form of an integral in each interval where du/ dr does 
not change sign, 

rr ru oj<£ 

r -r a = dt = {sign(dw/dr)} / (5) 

J r« } ^ p{V{u{ti))-V{x)) 
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In such conditions, the motion of a particle in a well is a periodic motion 
characterized by the function T\ 



71(^(0)) = f dX (6) 

7 « + (o) y/2{V{u+(0)) - V{x)) 

(For negative u(0) the integration limits shall be reversed). Note that T n 
may be expressed in terms of T x : 

T n (u + (0)) = (2[^] - 1)^(^(0)) + 2[^]T 1 («-(0)), (7) 

where [q] means the integer part of q, and u~(0) is defined from V(u + (0)) = 
V^(«-(0)) = E. 

For a given potential, the equation T n = R may be explicitly written and 
the existence of solutions to Problem (3 - 4) may be explicitly investigated. 

We are not going to give further details of the analysis in this simple case 
and turn out to the higher dimensional (D > 1) problem, i.e. motion with 
friction. In this situation, there is another exactly solvable problem: the 
motion in a quadratic potential. 

2.2 Motion in a quadratic potential (The linear equa- 
tion) 

We consider a quadratic potential V(u) = (l/2)Aw 2 . The equation of motion 
(3) takes the form 

d 2 u D — 1 du 

\u — . (8) 



dr 2 r dr 

d 



The solution of this Eq. with initial condition 4^(0) = is expressed as 



r l-n/ 2j 

D/2-i\(V^ r ), where J is the Bessel function [8]. It is important to 
note that the main properties of the solution may be understood simply from 
the invariance properties of Eq. (8). 

LEMMA: T n does not depend on u(0) and is proportional to A -1 / 2 . 

PROOF: The Eq. is invariant under a change in the scale of u, and also 
under the transformation 

r — > C r r, A — > C\\, 
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where C r = C x . 

According to this Lemma, the function T n (u(0)) takes a fixed value that 
depends only on A and n. Varying appropiately the parameter A (the po- 
tential), one may fulfil the requirement T n = R. The corresponding set of 
parameters , {A„}, define the eigenvalues of the linear problem. 

2.3 Some useful results 

In this Subsection, we derive a few general results following from the analogy 
with Mechanics and classify the potentials to be studied. 

In the presence of dissipation, the rate of change of the energy is written 

as 



i.e. u{r) is damped, as mentioned above. It means that E(u(0)) = V(u(0)) > 
E(0) > V(0) = (we have supposed that / is integrable, so that V(u) may 
be defined as / " f(x)dx). Then, we arrive at the following 

THEOREM (A necessary condition): If u(r) is a solution to (3 - 4) and 
/ is integrable, then u(0) is such that 



The last condition on the sign of f(u(0)) means that the particle shall 
be pushed towards the origin at the initial position, u(0). Otherwise, it will 
never move to the origing passing through u(0) because of the energy losses. 

More sophisticated versions of this Theorem will be formulated below 
when studying potentials with barriers. 

A second important result concerns the retardation effect of friction. Let 
us suppose that the particle moves from u a to Uf,. The time it spends in this 
motion may be written as 



dE/dr = ^-({\/2){du/drf + V(u) 



D - 1 



r 



(du/drf < 0, (9) 




f(x) dx > 0, sign(/(ii(0))) = sign(u(0)). 



/ 

J u 



dx 



r b -r a = 



^{dujdtaf + 2V(u a ) - 2V(x) - 2(D - 1) J ^{dx/drY ' 

(10) 
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Of course, this is not a closed expression because the derivative in the 
time interval (r a , r*,) enters the r.h.s. of it. However, it is evident that 



f u b ax , , 

n-r a > i 11 

Jun. . (An IA+ ^2 — I — OT/^, ^ _ OM(rr\ 



^f(du a /dt a y + 2V(u a )-2V(xy 
i.e. 

LEMMA: The time interval — r a is greater than the time the particle 
spends to move from u a to Ub without friction. 

Finally, let us classify the potentials according to their properties in the 
neighborhood of u — 0. In the present paper, we will study four classes of 
potentials having different behaviours in the vicinity of this point (Fig. 3): 

a) The wells are defined as concave potentials around u — 0. 

b) The hills are convex around u — 0. Of course, at large \u\, V(u) shall be 

positive (the necessary condition). 

c) and d) correspond to singular potentials. 

We will study below each class of potentials separately. 

3 Wells around u = 

A well is defined as a region with only one local extremum, the minimum at 
u — 0. In this Section, we study some examples of potentials having a well 
around u — 0. 

3.1 Potentials, quadratic in u = and \u\ — > oo 

Let V(u) be a potential such that 

V\ u ^o « ^A(0)« 2 , 
V\ H ^ « \K™)u\ (12) 



additionally, we will asume that the only zero of / is at u — 0. Then, we 
have the following 
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THEOREM: If A(0) < Ai and A(oo) > X k , then Problem (3 - 4) has at 
least 2k + 1 solutions. 

This Theorem was obtained in [5] . We will give a detailed proof of it by 
means of our method as an illustration. For the incomming Theorems, the 
proof will be shortened. 

The statement is that the function T n vs u(0) has the form depicted in 
Fig. 4 for 1 < n < k, i.e. for each T n there are two solutions. 

Indeed, the very small amplitude motion is governed by the u — > 
asymptotics of V . T n depends very smoothly on u(0) in this region and 
T n > T\ > R. The latter inequality comes from A(0) < Ai. 

On the other hand, the large amplitude motion is governed by the \u\ — > 
oo asymptotics and, according to the inequality A(oo) > A&, we have T n < 
T k < R. The point to clarify is why T n for large u(0) is not affected by the 
small-w behaviour of V. 

The answer is that, when -u(O) is large, the time the particle spends to 
move in the small-w region is negligible. This result comes from the scale 
invariance of the quadratic potential as shown in Fig. 5. Shadowed areas 
correspond to motion in the region \u\ < u a . It is seen that when |w(0)| — > oo 
the time spent in this motion shrinks to zero. It means that one can deform 
V(u) at low \u\ without changing significantly T n . 

Thus, Problem (3 - 4) has 2k nontrivial solutions plus the trivial u = 0. 

3.2 Potentials with barriers 

In the previous Subsection, we assumed continuity of T n vs u(0). However, 
continuity is broken when V has local extrema, others than u = 0. The point 
is that, as may be seen from Eq. 5 and the retardation Lemma of Section 
2.3, the time the particle spends to move out of a local maximum tends to 
infinity when u(0) approaches the position of the maximum. 

Then, let us first suppose that / has a unique second zero at a point 
a > 0. The following Theorem may be formulated 

THEOREM. If A(0) > A fc and VX (0) ^_oo > V(a), then Problem (3 - 4) 
has at least k solutions with < u(0) < a. 

To prove it, we draw again the function T n (u(0)), with 1 < n < k and 
positive u(0). At low u(0), T n < T k < R. When u(0) — > a from below, 
T n — ■> oo. The condition on V{— oo) guarantees that the motion is oscillatory 
around w(0) and the particle does not escape to — oo. 
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Note that it is difficult to draw the dependence T n vs u(0) for negative 
u(0) without a knowledge of the potential. The following Theorem, contained 
in Ref. [5], states that for asymptotically quadratic potentials one can say 
much more. 

THEOREM. If / has positive zeroes, the first of which is at u = a, and 
A(0), A(oo) > Afc, then Problem (3 - 4) has at least 4k — 1 solutions. 

We have drawn in Fig. 6 the potential and the functions Ti, T n , 1 < n < 
k. The points b + and 6_ are defined in the monotone regions. They satisfy 
V(a) = V(b + ) = V{b). Dashed lines means that the curves are conditionally 
drawn, while shadowed intervals of w(0) mean physically impossible initial 
conditions. 

The dependence Ti on -u(O) when < w(0) < a is the same as in the 
previous Theorem. For very large positive w(0), Ti is determined by A(oo), 
i.e. TiI^q,^ < Tj, < R. On the other hand, because of energy losses, if 
the particle starts from b + it will not reach the origin. By continuity, there 
exists c\ > b + such that the particle arrives at a with zero velocity. This 
corresponds to an infinite T\. When u > c\ the particle reaches the origin 
and the dependence Ti(w(0)) is shown. Note that we can not say anything 
about Ti for negative -u(O). Thus, the equation Ti = R will have, at least, 
two solutions. 

Analog reasonings are used in the analysis of T n , 1 < n < k. C n is now 
defined such that when -u(O) > c n the origin is reached n times. Note that 
T n (c„) = oo and also that c\ = C2 < c 3 = C4 < c 5 . . .. On the l.h.s. of the 
origin, we can define the points e n < d < b_. d is such that when the particle 
arrives to a it do so with zero velocity, while e n is such that for -u(O) < e n , 
the particle reaches the origin n times. Note that €2 = e 3 > = e§ > e§ . . .. 
In other words, for each n there are 4 solutions. This proves the Theorem. 

Notice that, unlike papers [1 - 5], we are able to indicate forbidden regions 
for w(0). This is a generalization of the necessary condition of Section 2.3. 

3.3 The potentials V = g\u\ p 

Let us now consider the potentials V = g\u\P, with g > 0, f3 > 1. We shall 
first prove that, whatever f3 be, u(r) will have the form drawn in Fig. 2. 
After that, we will use scale- invariance properties of the equation of motion 
to obtain the dependence T n vs u(0). Let us prove the following general 
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LEMMA. In a potential well, u(r) is an oscillating function of decaying 
amplitude. 

PROOF. It is evident that the particle will reach the origin whatever the 
initial position be. It can not stop in an intermediate point where the force 
is not zero. Thus, the question is how long it takes to reach the origin and 
what is the final velocity. If this time and the velocity are finite, we can 
repeat the argument to conclude that u(r) will have infinite zeroes. 

Let r a be an intermediate time such that ^(r a ) > 0. Due to the par- 
bound for the time to 

du i 
dr 1 



dr 

ticular form of the friction, we can obtain an upper 
reach the origin starting from u(r a ), r^, and a lower bound for 4^(r&) , if we 



neglect the potential for r > r a and solve the problem: 

d 2 u D — 1 du 



dr 2 ~^ r dr ' ^ ^ 

u(r a )=u a , - r (r a )=v a , (14) 



dr 

which has the following solution 

^(r) = vMr) - 1 , (15) 
«(r) = u a + v a ^{rl- D -r 2 - D } . (16) 

It means that u — will be reached in a finite with a finite velocity 
and u(r) will have infinite zeroes. 

Thus, let us now turn out to the dependence T n vs u(0) in the potentials 
V = g\u\P. The equation of motion takes the form 

_ = -^(sign u)\uf ^ • (17) 

The properties of T n following from the scale invariance of the equation 
are given in the next Lemma: 

LEMMA. For fixed g, T n ~ ^(O)! 1 ^/ 2 , while for fixed u(0), T n ~ g- 1 ' 2 . 

Thus, for every n, the equation T n = R will have two solutions, and we 
arrive to the following 
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THEOREM. Problem (3 - 4) with V = g\u\ p , g > 0, (3 > 1 has infinite 
solutions [3]. 

One can now combine these with previous results. In quality of example, 
let us formulate the following 

THEOREM. Let A(0) < X k and VL ^ ~ \u\ p , with (5 > 2, then solu- 
tions to Problem (3 - 4) with any n > k zeroes exist. 

The curve T n vs u(0) for n > k may be easily drawn in this case. Note 
that the dependence of Tk|| u (o)|^oo on the low-u(O) properties of V is, for 
(5 > 2, weaker than in the quadratic potential. 

4 Hills around it = 

We now study the motion in a potential like that one shown in Fig. 3 b. 
For simplicity, we assume that V is quadratic near zero (A(0) < 0) and also 
quadratic at large values of u. No additional zeroes of / exist. Then one can 
formulate the following 

THEOREM. If A(0) < and A(oo) > A fc , then Problem (3 - 4) has 2k + 1 
solutions. 

We have drawn in Fig. 7 the curve T n vs u(0) for 1 < n < k. The 
large-w(O) behaviour of it is evident. The points b + and 6_ are the zeroes of 
V. The points c n and e n are defined as in the previous Section, i.e. starting 
from the right of c n (the left of e n ) the particle may reach the origin n times. 
Note that would it start from c n (e n ), then it would arrive to u — with 
zero velocity, i.e. T n (c„) = T n (e„) = oo. Note also that b + < c\ < c 2 • • •, 
6_ > ei > e 2 > • • •. Thus, for each n there are two solutions and the Theorem 
is proved. 

Other potentials could be analysed, but we think that to show the ad- 
vantages of the method the given example is enough. 

5 Singular Potentials 

The main property of the singular potentials, Figs. 3 c) and d), is that the 
force, — dV/du, at u = is ill-defined. So, the motion an interval of time after 
the particle reaches the origin is not well defined, and we can only analyse 
the existence of positive solutions to (3-4). 
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An example of a potential like 3 c) is V = g\u\", with g > and < /3 < 1. 
Let us stress that the upper bound for and the dependence Ti ~ | w(0) 1 1_/3 ^ 2 , 
obtained in the Lemmas of Section 3.3, are valid, so that the equation Ti = R 
has always a solution in this case. 

The same analysis holds for the potential V = — g|w|~ /3 , with g, f3 > 0. 
This is a potential of the form 3 d). Scale invariance in this case leads to 
Ti ~ |w(0)| 1+/3 / 2 , so that the equation Ti = R will always have a solution 
also. 

We can now combine possibilities to obtain interesting situations. Let, 
for example, the potential V be quadratic at the origin with A(0) > 0, while 
at long distances V ~ Vo — g\u\~P. No zeroes of / exist, except the trivial at 
u = 0. Then, we obtain the following 

THEOREM. If A(0) > A fe , then Problem (3 - 4) has at least 2k + 1 
solutions. 

The proof is trivial. 

6 Concluding Remarks 

In the present paper, we used the analogy of Eq. (3) with the second Newton's 
law in order to obtain existence theorems in Problem (3-4). Appart from 
reproducing existing results, we give new examples of potentials (of /) in 
which it is relatively easy to analyse the existence of solutions. 

We think that the given examples show that the method is general enough 
to provide a first insight to the problem for any reasonable function /. After 
that, we may go further on in two ways: i) Use more rigurous methods to 
complete the proof and/or ii) Obtain numerical solutions to the equation. 
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Figure Captions 



Fig. 1 

a) The analogy with Mechanics. 

b) A positive solution to (3 - 4) corresponding to the situation depicted 
in a). 

Fig. 2. A generic damped oscillating function u(r) describing the motion 
of a particle in V. 



Fig. 3. Different possibilities for the neighborhood of u = 0. 

a) Well 

b) Hill 

c) Finite, but singular 

d) Infinite, singular potential 

Fig. 4. Dependence T n vs w(0) for the potential considered in Section 

3.1. 



Fig. 5. A consequence of the scale invariance of the quadratic potential. 
The shadowed areas correspond to motion in the region \u\ < u a . When 
|w(0)| — > oo, the time spent in this motion shrinks to zero. 

Fig. 6. A potential with barriers and the corresponding Ti(u(0)), T n (w(0)), 
1 < n < k. 

Fig. 7. The curves T n vs w(0), 1 < n < k, for the potential of Section 4. 
Notations are the same as in Fig. 6. 
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